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(57) Abstract 

A carrier tracking loop that uses a frequency locked loop (FIX) and incorporates synchronized matched filters to minimize complexity 
of the system. The synchronized matched fitters eliminate out-of-band interference while minimizing computation. The earner tracking 
loop has an improved error signal generation to extend a frequency range and by employing error signal [normalization, which extends the 
dynamic range of the amplitude of the input signal and allows the frequency to "flywheel" when signal fading occurs. 
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CARRIER TRACKING LOOP FOR OPSK DEMODULATOR 
Related Applications 
U.S. Patent Application entitled "Mobile 
Terminal Apparatus and Method for a Satellite 
Communication System" filed by Albert: J. Fasulo II, 
5 Denise M. Cammarata, Keith W. Janson, Samuel S. Anderson, 
Raymond R« Cooper, and Roy Stehlik, on the same date as 
the instant application which is a continuation-in-part of 
U.S. Patent Application Serial No. 08/242,400 filed 
May 13, 1994, is incorporated by reference. 
10 U.S. Patent Application entitled "Carrier 

Acquisition Technique For Mobile Radio QPSK Demodulator" 
filed by Brian W. Kroeger, Joseph B. Bronder, and 
Jeffrey S. Baird on the same date as the instant 
application. 

15 U.S. Patent Application entitled "Locally 

Coherent QPSK Detection With Differential Decoding For A 
Fading Channel" filed by Brian W. Kroeger, Joseph B. 
Bronder, and Jeffrey S. Baird, on the same date as the 
instant application, is incorporated by reference. 

20. U.S. Patent Application entitled "A Digital 

Quadriphase-Shift Keying Modulator" filed by Brian W. 
Kroeger and Roy Stehlik, on the same date as the instant 
application, is incorporated by reference. 
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U.S. Patent Application entitled "Numerically 
Controlled Oscillator With Coupled Exponential Outputs 
Using Recursion Technique" filed by Brian W. Kroeger and 
Jeffrey S. Baird, on the same date as the instant 
5 application, is incorporated by reference. 

U,S. Patent Application entitled "Symbol 
Synchronizer Using Modified Early/ Punctual /Late Gate 
Technique" by Brian W. Kroeger, Joseph B. Bronder, Tod A. 
Oblak and Jeffrey S. Baird, on the same date as the 
10 instant application, is incorporated by reference. 

Background of th «» Tnvpntion 
The present invention relates to a mobile 
terminal for a communication system, and more particularly 
to a method and apparatus for digital receiving in a 
15 communication system. 

While the invention is subject to a wide range 
of applications, it is especially suited for use with a 
satellite communication system, and will be particularly 
described in that connection. 
20 Description of Related Art 

In a typical satellite communication system, 
outgoing RF signals transmitted from a mobile terminal 
unit are received directly by the satellite. The 
satellite in turn retransmits the RF signals to a ground 
25 station that is connected by wire to a public switched 
telephone network (PSTN) , which in turn routes the 
outgoing signals to either a conventional telephone or to 
another mobile terminal unit of a satellite or cellular 
network. Incoming signals from a conventional wired 
30 telephone are conducted from the PSTN to the satellite 

ground station, which in turn transmits -RF signals to the 
satellite for retransmission to the mobile terminal unit. 
Thus, communication can be between two mobile terminal 
units or between a mobile terminal unit and a conventional 
35 telephone connected to a PSTN, or between a satellite 

mobile terminal unit and a cellular mobile terminal unit, 
for example. In each of the aforesaid conditions (except 
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mobile terminal to mobile terminal ) , the communication is 
routed through a PSTN. 

It has been proposed, to provide a satellite 
communications network that utilizes geosynchronous 
5 satellites each of which is capable of covering an area 
corresponding to a substantial portion of the North 
American continent, so that one satellite with 
approximately six beams will cover the entire continent 
from Alaska to Mexico. The satellites for such a network 

10 will be approximately 22,600 miles above the equator and 
will be designed to operate in the L-band of RF 
frequencies. For example, the frequency of the signal 
being transmitted to the satellite will be between 
1626 MHz to 1660 MHz, and the frequency of the signal 

15 received from the satellite will be between 1525 MHz to 

1559 MHz. Energy travelling this great distance undergoes 
huge attenuation such that the power flux density incident 
at the antenna of the mobile unit is approximately 10~ 14 
watts per square meter. This grossly attenuated signal is 

20 further degraded by noise, and other satellite channel 

impairments such as Rician fading due to multipath-f ading. 

Mobile terminal units capable of receiving 
signals of this order of magnitude are either limited to 
paging signals or require extremely sophisticated 

25 hardware , which is bulky and heavy and does not lend 
itself to mobility. Also, the limited functions and 
expense of such terminal equipment would prevent its 
acceptance among a large segment of potential users. 

In light of the foregoing, there is a need for a 

30 mobile terminal unit that is not only capable of reliably 
receiving the attenuated signals, subject to Rician 
fading, of a high altitude satellite for voice, facsimile, 
and data communication, but also compact, lightweight, and 
relatively inexpensive to manufacture. 

35 gflffTHHTY the Invention 

The present invention overcomes the problems and 
disadvantages of the prior art by using a frequency locked 
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loop (FLL) in place of a conventional phase locked loop 
(FLL) and by incorporating synchronized matched filters to 
minimize complexity of the system. The synchronized 
matched filters eliminate out-of-band interference while 
5 minimizing computation. The filter outputs are also used 
by other parts of the system. In addition, the present 
invention overcomes the problems and disadvantages of the 
prior art by including an improved error signal generation 
to extend a frequency range and by employing error signal 
10 normalization, which extends the dynamic range of the 

amplitude of the input signal and allows the frequency to 
"flywheel" when signal fading occurs. 

In accordance with the purpose of the invention, 
as embodied and broadly described herein, the invention is 
15 a demodulator for a mobile telecommunications system, 

including a carrier tracking loop receiving a series of 
symbol samples, each having an inphase signal and a 
quadrature signal, comprising: a complex multiplier 
receiving the sequential received samples? a matched 
20 filter pair, connected to the complex multiplier and 
outputting a complex filtered signal; means for 
determining a phase difference between two sequential 
received samples; means for determining a normalized angle 
in accordance with the phase difference; means, connected 
25 to the normalized angle determining means, for determining 
an increment phasor from the normalized angle; and 
means, connected to the increment phasor determining means 
and to the complex multiplier, for generating a complex 
exponential value that is input to the complex multiplier, 
30 wherein the complex multiplier includes: means for 

inputting the complex exponential value to the complex 
multiplier, and means for performing a complex 
multiplication operation on the complex exponential value 
and a current pair of the complex inphase /quadrature 
35 signals to translate the complex ihphase/quadrature 
signals to approximately baseband frequency. 



WO 96/12367 



PCT/US95/11950 



5 

Objects and advantages of the invention will be 
set forth in part in the description which follows and in 
part will be obvious from the description or may be 
learned by practice of the invention. The objects and 
advantages of the invention will be realized and attained 
by means of the elements and combinations particularly 
pointed out in the appended claims . 

Qyief Description of tfre Drawings 

The accompanying drawings , which are 
incorporated in and constitute a part of this 
specification, illustrate several embodiments of the 
invention and, together with the description , serve to 
explain the principles of the invention. 

Fig. 1 is a schematic block diagram of a mobile 
terminal apparatus ; 

Fig. 2 is a schematic block diagram of a 
demodulator of the system of Fig. 1; 

Fig. 3 is a schematic block diagram of a carrier 
tracking loop for the demodulator of Fig. 2; 

Fig. 4 is a chart showing an encoding scheme for 
a locally coherent QPSK system; 

Fig. 5 is a schematic block diagram of a 
recursive numerically controlled oscillator of the carrier 
tracking loop of Fig. 3? 

Fig. 6 is a graph showing a relationship between 
a value of a phase increment 0 input to a numerically 
controlled oscillator of Fig. 3 and a cosine value output 
from the oscillator; 

Fig. 7 is a graph showing a performance of the 
numerically controlled oscillator of Fig. 5; 

Fig. 8 is a schematic block diagram of a matched 
filter pair and a symbol synchronizer of Fig. 3; 

Fig. 9 is a schematic block diagram showing more 
detail of the matched filters of Fig. 8; 

Fig. 10 is a graph of the output of a normalizer 

of Fig* 8; 
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Fig. 11 is a graph of output of a first 

integrator of Fig. 8? 

Fig. 12 is a graph of output of a second 

integrator of Fig. 8; 

Fig. 13 is a graph of an interpolation control 

signal of Fig. 8; 

Fig. 14 shows a exemplary frequency spectrum of 
an incoming signal to the carrier tracking loop of Fig. 3 
before it is shifted to near baseband frequency; 

Fig. 15 shows a frequency spectrum of the signal 
of Fig. 14 after it is shifted to near baseband frequency; 

Fig. 16 shows an example input to the matched 
filter pair and symbol synchronizer of Fig. 3; 

Fig. 17 shows an example of the tinning of the 
matched filter leading the timing of the incoming signal; 

Fig. 18 shows an example of the timing of the 
matched filter lagging the timing of the incoming symbol; 
and 

Fig. 19 shows examples of the timing of filter 
coefficients in the early, punctual, and late filters. 
no^ilgd Description of Preferred Embodiment 

Reference will now be made in detail to the 
preferred embodiment of the invention, examples of which 
are illustrated in the • accompanying drawings. Wherever 
possible, the same reference numbers will be used 
throughout the drawings to refer to the same or like 
parts. 

Fig. 1 shows a mobile terminal unit comprising a 
transceiver 12, an antenna 14, and one or more user 
interfaces 16. The transceiver 12 has an RF section 18 
and a processor board 20. 

The RF section has a receiving portion that 
includes a low noise amplifier 22, a down converter 26, 
and an analog to digital converter (A/D) 28. The RF 
section also has a transmitter portion that includes a 
high power amplifier (HPA) 30, an upconverter 32, signal 
line 62, and digital-to-analog converter (D/A) 34. The 
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processing is performed by the DSP 38 and 39 as controlled 
by the control processor 40. Received signals are 
demodulated, in other words translated, from their 
waveform representation into their digital data 
5 representation, decoded and packaged for voice, fax or 
data processing at the DSP 39. 

The processor board comprises a signal processor 
including a digitally implemented demodulator function, a 
digitally implemented modulator function, a digitally 
10 implemented receive framing function coupled to the 

demodulator, a digitally implemented transmit framing 
function coupled to the modulator function, and a 
. digitally implemented encoding function and decoding 
function coupled to the transmit and receive function, 
15 respectively. 

As the QPSK signal is received, several types of 
synchronization must be performed on the signal. First, 
the system must determine the frequency of the "carrier 
wave" that is being phase modulated. Second, the- system 
20 must determine where each symbol (represented by multiple 
bits) begins, i.e., "symbol synchronization." 

Fig. 2 is a schematic block diagram of a 
demodulator 264, which is implemented as instructions 
executed in DSP 38 of the system of Fig. 1. In Fig. 2, an 
25 analog signal 205 is received and converted to a digital 
signal 206 by the analog-to-digital (A/D) converter 28. 
Signal 206 is input to Hilbert Transformer 278, which 
yields a digital 10-bit in-phase signal (I) 207 and a 
digital 10-bit quadrature signal (Q) 208. Signals (I) 207 
30 and (Q) 208 are input to carrier tracking loop 220, which 
includes a complex multiplier 230, matched filters and 
symbol synchronizer 240, a symbol synchronizer 
estimator 250, a carrier tracking element 260, and a 
carrier frequency estimator 270. The output from the 
35 carrier tracking loop 220 is input to a QPSK detector 280. 

The elements of demodulator 264 preferably are 
embodied in instructions performed by processor 38 of 
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antenna 14 is connected to the RF section through a 

diplexer 36. 

The processor board has a digital signal 
processor 37 including a modem processor 38 and a 
coder /decoder and fax processor 39, a control 
processor 40, a universal asynchronous 

receiver/transmitter (UART) 42, an A/ D converter 44, a D/A 
converter 46, and a fax interface unit 48. The user 
interfaces of the terminal unit include an output 52 for a 
handset, an output 54 for a facsimile machine and an 
output 56 for connection to a personal computer. The 
apparatus may also include a cellular transceiver board. 

The transmit signal begins in the processor 
board 20 as digital data, which is then converted into 
digitized, modulated waveform samples. This digitized 
waveform, which is created by the processor 37, is 
conducted to the digital-to-analog (D/A) converter 34 
before the upconverter 32 of the RF section. The 
converter 34 converts the digitized waveform into both an 
analog I (in phase) and analog Q (quadrature) waveform. 
The analog I and Q waveforms are then sent to the mixer or 
upconverter modulator 32 where they are upconverted to an 
L-Band, and then summed to form a low level L-Band 
transmit signal on line 62. This low level L-Band 
transmit signal is amplified by the high power 
amplifier 30. The resulting high level L-Band transmit 
signal then passes through transmit side of the 
diplexer 36 and is steered towards the MT antenna 14 from 
which it travels to the satellite 60. 

In operation, the DSP 38 reads the received 
signal waveform samples from the A/D converter 28, and 
processes these samples as they arrive. The L-Band 
receive signal is mixed down to the low IF frequency to 
aid in signal filtering and to reduce the A/D sample rate 
so that the digital signal processing can keep up wxth the 
sample arrivals. Once the receive signal has been sampled 
and loaded into the DSP 38, all subsequent signal 
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Fig* 1, which executes a software program stored in a 
memory, such as a ROM (not shown). Processor 38 
preferably is a DSP 320C51 fixed point processor 
manufactured by Texas Instruments Corporation. However, 
any processor capable of performing the process described 
could be used to implement the invention. The bit-widths 
of the signals in the carrier tracking loop 220 are chosen 
to provide sufficient dynamic range for the particular 
application by minimizing the effects of quantization 
noise while preventing numerical overflow. 

Certain elements of the demodulator 264, such as 
oscillator 314, symbol synchronizer 240, or symbol synch 
estimator 250 may be implemented in hardware. A person of 
ordinary skill in the art understands that Fig. 2 can 
represent either a hardware or a software implementation 
of the present invention, and that the invention is not 
limited to either a hardware or a software implementation. 

Although carrier tracking loop 220 is shown to 
include oscillator 314, symbol synchronizer 240, and 
symbol synch estimator 250, it should be understood that 
carrier tracking loop 220 may also be implemented using 
other, conventionally known elements in place of any or 
all of at least elements 314, 24 0, and 250. In addition, 
oscillator 314 may used in other systems requiring an 
oscillator and is not limited to usage in the system shown 
herein. Filter element 240 and/or symbol synch 
estimator 250 also may be used in other systems and are 
not limited to usage in the system shown herein. 

In accordance with the present invention, the 
QPSK detector uses a technique for locally coherent 
detection (demodulation) of a digital QPSK modulated 
signal that performs especially well in a M fading 
environment • ** The parameters of the technique can be 
optimized to suit the particular fading characteristics 
and statistics in order to optimize performance in a 
particular known fading environment* The locally coherent 
QPSK detector 280 of the present invention produces a 
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multi-bit "soft-decision" output that is then input to a 
forward error correction (FEC) element (not shown) . ■ 

Fig. 4 is a chart showing an encoding scheme for 
a locally coherent QSPK system. Fig. 4 shows how bit ^ 
pairs are encoded. The encoder output is a complex number 
whose phase is 45 degrees plus an integer multiple of 
90 degrees. The encoder input determines the difference 
between the phases of two consecutive encoder _ 
Thus, if, for example, the encoder input is 1,1, then the 
> difference between the phase of the current encoder output 
and the previous encoder output is always 180 degrees 
independent of the actual value of the previous output. 
Similarly, the phase difference between current and 
previous encoder outputs is 90 degrees when the encoder 
5 input is 0,1. Thus, the phase differences of the encoder 
. outputs contain the information necessary to recover the 
bit pairs that make up the encoder inputs. 
I. A Carrier Tracking Loop in a 
Mobile Communic t^n System 
20 Fig . 3 is a schematic block diagram of the 

carrier tracking loop 220 for the demodulator 264 of 
Fig. 2. The carrier tracking loop 220 includes the 
complex multiplier 230, the matched symbol 
synchronizer 240, the symbol synchronization 
25 estimator 250, a fourth power element 302, the carrier 
tracking element 260, and a numerically controlled 
oscillator (HCO) 314. Carrier tracking element 260 
includes a frequency discriminator 304, a frequency error 
estimator 306, a gain element 308, a lossy digital 
30 integrator 310, and a perfect integrator 312. In 

implementation an initial value of integrator 312 is set 
in accordance with a signal input from the carrier 
frequency estijaator 270 of Fig. 2, where the input 
represents an initial frequency estimate from 

35 estimator 270. 

The purpose of the carrier tracking loop is to 

minimize the carrier frequency offset error after the QPSK 
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modulated signal has been translated from a radio 
frequency or an intermediate frequency to near baseband 
frequency. Pig. 14 shows a frequency spectrum of an 
incoming signal before shifting. Fig. 15 shows a 
5 frequency spectrum of the spectrum of Fig. 14 after it is 
shifted by the complex multiplier 230. 

The carrier tracking loop 220 is a closed-loop 
servo mechanism which , through negative feedback, attempts 
to keep the frequency of the NCO 314 equal to the complex 

10 conjugate of the incoming offset-baseband carrier 

frequency. Loop 220 is a frequency locked loop (FLL) . 
The digital 10-bit in-phase signal I 207 and the digital 
10-bit quadrature signal Q 208 are input to the complex 
multiplier 230. Output signals 281 and 282 from the 

15 matched filter pair 240 have some arbitrary phase 

rotation, some frequency offset, and a frequency tracking 
error . 

The carrier tracking loop 220 is designed with 
parameters a and 13 such that residual frequency - 

20 perturbations are sufficiently small to minimize adverse 
effects on the subsequent detection by QPSK detector 280. 
Parameters a and R are chosen to yield a mini m u m frequency 
tracking response time (a fast reacquisition time) with 
minimum M ringing" and a low loop noise bandwidth. The 

25 arbitrary phase offset (stationary or slowly changing) is 
removed by the subsequent QPSK detector 2 80. 

The offset baseband signal 207 and 208 input to 
the complex multiplier 230 of the carrier tracking 
loop 220 can be represented as follows: 

30 (1) 

A 16-bit sine signal 320 and a 16-bit cosine 
signal 322 from the NCO 314 are also input to the complex 
multiplier 230. The complex multiplier 230 includes four 
multipliers 324 , an adder 326 f and a subtracter 328. The 
35 complex multiplier 230 outputs a digital 16-bit in-phase 
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signal I 330 and a digital 16-bit quadrature signal Q 332. 
After mixing (i.e., multiplication by the complex 
multiplier 230) the complex number represented by I,Q 
signals 330 and 332 are input to the matched filter pair 
5 and symbol synchronizer 240. 

In the matched filter pair and symbol 
synchronizer 240, the frequency corrected baseband signal 
is filtered at an optimum symbol timing. The output of 
the matched filter 240 is also required for operation of 
10 the QPSK detector 280 of Fig. 2 and are output thereto. 
By placing the matched filter 240 within the carrier 
tracking loop 220, out of band signals are filtered out, 
allowing the loop to track with reduced interference. The 
output symbol samples from the matched filter 240 have 
15 arbitrary phase offset and a small residual frequency 
perturbation (frequency tracking error) that are not 
corrected by the carrier tracking loop. Such an error is 
due, e.g., to acceleration or drift from oscillator 314. 

The NCO 314 generates a complex exponential 
20 (sine and cosine signals 320 and 322), which is as close 
as possible to the complex conjugate of the incoming 
signal 207 and 208. 

(2) 

In closed loop operation, the frequency 
25 perturbation (Af ) term of equation 2 is a zero mean random 
process. The input to the matched filter element 24 0 is 
then: (3) 

The input sample rate with index m is greater 
30 than the symbol rate with index n by a preferred factor 

equal to the number of complex input samples representing 
each symbol. In a preferred embodiment, eight samples are 
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taken for each symbol. The symbol synchronization output 
samples S^fn) of the symbol synchronizer 24 0 preferably is 
represented by: 

(4) 

where d is 0,1,2,3 dependent upon the pair of information 
bits representing the nth QPSK symbol, T s is the symbol 
period, and ^ n is the arbitrary phase offset • The digital 
8-bit in-phase signal I 281 and the digital 8-bit 
quadrature signal Q 282 are output from the matched filter 
pair 240 to QPSK detector 280 of Fig. 2. The output rate 
of the symbol synchronizer 24 0 is 1/8 the rate of the 
incoming sampled symbols. 

Signals I,Q 281 and 282 also are input to fourth 
power element 302. Fourth power element 302 removes the 
modulation in the signal whose phase constellation is 
spaced at tt/2. Since the fourth power multiplier the 
exponent by four, then 4-7T/2 is 2tt, which is the modulus 
of the phase. Thus, the fourth power element 302 
eliminates modulation in the signals 281 and 282: 

(5) 

[5^(i3)] 4 = a 4 •e : ' (8 ** A^ ■ flT »* 44 " , +noise4(n) 

The fourth power signal is then input to the 
frequency discriminator 304 of the carrier tracking 
element 260, which computes a complex difference phasor 
between successive samples. Frequency discriminator 304 
inputs a current sample and includes two one-period delay 
elements 340, which yield a previous sample. Frequency 
discriminator 304 performs a complex multiply operation on 
the current symbol and the previous symbol using four 
multipliers 342, the adder 344, and the subtractor 346. 
The output of frequency discriminator 304 can be 
approximately represented (approximation assumes 
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10 



15 



20 



30 



successive 
by: 



samples have equal magnitude and phase offset) 



(6) 



[S^n) ]«■< [S^a-D ]*)' -a* a -e^+ fT '+noise8(n) 



Frequency discriminator 304 outputs X signal 350 and Y 
signal 352 to the frequency error estimator 306 of the 
carrier tracking element 260. 

The frequency error estimator 306 determines an 
instantaneous frequency error 8. The frequency error 
estimator 306 includes two absolute value elements 354 and 
356, a subtracter 358, two adders 360 and 362, a sign 
determiner 364, a multiplier 366, and a normalizer 367. 
In other implementations of the invention, the .frequency 
locked loop of carrier tracking loop 220 derives an error 
signal 6 from the imaginary element (I) of the previous 
expression, or uses y/abs(x) to -normalize" the dynamic 
range of the output from carrier tracking element 260. 
Both of these techniques, however-, limit the frequency 
range of the phase discriminator, since the error signal 
starts to diminish as the magnitude of the frequency error 
exceeds 1/16 of the symbol rate. In another 
implementation, the frequency range is extended {and 
response time improved) by computing the arctangent of 
(x+jy) over 2n radians. This embodiment extends the 
frequency range of the frequency locked loop to -1/8 of 
25 the symbol rate. However, the 2rr arctangent is 

computationally complex, and responds to noise when the 
signal fades. Thus, various ones of the described 
implementations are used, depending on the operating 

environment • 

in Fig. 3, the frequency error estimator 306 
receives X signal 350, Y signal 352, and a constant NORM 
signal 355. The frequency error estimator 306 determines 
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an instantaneous frequency error in accordance with the 
following equation: 

(7) 

e _ y n +(\X n \-X n ) -sgn{y n ) 

\X4+\Y n \+NORM 

In the absence of noise, an error signal 9 n 368 
is a monotonically increasing function of the arctangent 
over +-rr radians. In the presence of noise , the output 8 n 
368 of the frequency error estimator 3 06 is actually 
represented by the same function convolved with the 
probability density function of the phase noise process. 
Furthermore, in the presence of noise, the error signal is 
suppressed as the magnitude of signals X and Y 350 and 352 
drop below the constant NORM 355. NORM 355 is preferably 
3dB above the noise floor (i.e., the mean of (iX[+|Yj) in 
the absence of a signal). 

Use of the NORM signal 355 allows the carrier 
tracking element 260 to "flywheel" through signal fades. 
The expression allows the input signal to operate over a 
very high dynamic range without changing the 
characteristics of the carrier tracking loop 220. The 
NORM signal acts to "suppress" output from the carrier 
tracking element 260 when signal fade occurs. This is in 
contrast to conventional systems, in which the loop gain 
increases with the signal magnitude. The value of NORM is 
approximately set to the magnitude of x+jy corresponding 
to the point where the signal approaches a fading 
condition and the loop response is to be suppressed. 

The NORM signal 355 preferably is determined as 
described in the following paragraph. If the noise 
variance of each symbol out of the matched filter is q 2 /8 
(e.g., 8 samples per symbol and q is set to the root mean 
squared (rms) thermal noise into the carrier tracking 
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loop 220), then the amplitude of the signal element 
(ignoring noise) of the x+jy signal is approximately: 



a-- 



The NORM signal 355 is set equal to this value 
at the appropriate level of Eb/N G where fading is 

determined to begin. 

The error signal 8 n 36.8 is then processed by 
the loop filter formed of, e.g., the constant gain a 
element 308, the lossy digital integrator 310 with gain A 
372, and the perfect integrator 312. The values of a and 
fl are chosen to establish the best compromise between 
response time, tracking accuracy, and ensuring closed-loop 
atability. The output of the perfect integrator 312 
controls the frequency of the NCO 314 as described below. 

in one implementation, a separate frequency 
acquisition algorithm is used to estimate the initial 
carrier offset error of the perfect integrator 312 and to 
set the initial frequency stored in the integrator 312, 
which provides the frequency control for the NCO 314. 
Because this frequency error can change over time due to, 
e.g., vehicle acceleration, the frequency offset error 
must be tracked and compensated. 

H. Tfro mi-ric ^W r-nntrolleH OttcjUator. 

Fig. 5 is a schematic block diagram of the 
recursive numerically controlled oscillator (NCO) 314 of 
the carrier tracking loop of Fig. 3. The NCO 314 uses a 
recursive method to generate complex exponential (sine and 
cosine) outputs 320 and 322. The described embodiment is 
a very efficient implementation of the invention and 
achieves high spectral purity of the synthesized 
exponential. Specifically, the phase noise is due only to 
the precision of the processor 38 that implements the 
computer instructions. Fig. 6 shows that as * increases, 
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a frequency of cosine signal 322 (and of sine signal 320) 
also increases • 

A sample plot of the spectrum output from the 
NCO 314 for a given <p is shown in Fig. 7. This plot shows 
5 a 1024 point FFT of a complex exponential generated by the 
NCO 314 at near the predicted worst case frequency. The 
exact frequency was chosen so that the side lobes were 
eliminated by uniformly windowing a whole number of cycles 
contained in the 1024 points. Integrating over the 1024 
10 frequency bins of the FFT, the total noise is about -120 
dB from the signal (60 dB or better than the predicted 
bound). 

As shown in Fig. 5, the NCO 314 includes a 
complex phasor approximator 502 , a complex multiplier 504 , 

15 and a recursive normalizer 506. The complex phasor 
approximator 502 inputs phase increment <p 39 0 from 
integrator 312 of Fig. 3 and outputs complex phasor 
increment 518 and 519. The complex phasor 
approximator 502 also includes a squaring element 510, a 

20 multiplier 512, a subtractor 514, and a delay element 516. 

The complex phasor approximator 5 02 converts the 
scalar phase increment signal 390 to a complex phasor 
increment 518 and 519. The scalar phase increment 390 is 
limited to an alias frequency that is conservatively less 

25 than that limited to about n/A radians per sample. During 
this conversion, processor 38 makes an approximation to 
the imaginary element (sine), then determines an 
appropriate real (cosine) element to optimize spectral 
purity. The following paragraph discusses how the 
.30 imaginary signal 519 and the real signal 518 are 
determined * 

For small values of <p (e.g., <p<l) , sin(^) * 
Thus, in Fig* 5, complex phase approximator 502 yields an 
unnormalized imaginary signal 519, the phase increment 390 
35 is passed through the delay element 516 (to compensate for 
the time used by the remainder of the complex phasor 
approximator 502), but is not otherwise altered. A person 



WO 96/12367 



PCTAUS95/11950 



18 



of ordinary skill in the art will understand that other 
delays (not shown) may occur in the carrier tracking 
loop 220 and specifically in the NCO 314 as needed. A key 
factor in producing a "pure" sine wave is to start wrth an 
unnormalized increment phasor that is as close to unity 
magnitude as is practical. The real element 518 is found 
by approximating a "cosine-like" function. It is 
desirable to obtain spectral purity, but not necessarily 
to obtain an exact resulting frequency. The frequency can 
be corrected by other means as described below. Thus, the 
real element 518 can be found such that the magnitude of 
the increment phasor is near unity. Instead of computing 
a cosine approximation, processor 38 expands a series to 
approximate the unity magnitude solution, resulting in the 
expression for real element signal 518: 



(9) 



a a +<t> 2 «i; 

where a is the real element 518 of the increment phasor. 
When this series is truncated, a positive bias exists. 
The error due to this bias can be reduced by adjusting the 
last coefficient in the truncated series. After 
truncation of all but two terms, a good approximation of 
the series for the real element 515 of the phasor 
increment is : 

a « 1-0. 62<p 2 

Therefore, the phasor increment approximation 518, 519 is 
defined as: (10) 

phasor_inc * a Us) ♦j-Jb(n) * i-o.62*£+j*,. 
where j is the square root of -1. The increment phasor of 
approximately unity magnitude is a function of the phase 
increment 390. 
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The complex multiplier 504 inputs the real 
signal 518 and the imaginary signal 519 and outputs real 
and imaginary signals 531 and 532- The complex 
multiplier 504 also inputs recursive feedback signals 550 
and 552 output from the recursive normalizer 506, The 
complex multiplier 504 includes four multipliers 520 , 522 , 
524 , and 526 f a subtractor 528 , and an adder 530- The 
complex multiplier 504 multiplies the phase increment 518 
and 519 and determines the previous complex exponential 
output . This is the basic recursive operation of the 
NCO 314. 

(11) 

c(n)+j-d(n) = [a(n-l) +j-2>-(u-l) 1 -phsLSorJnc 

where c(n) and d(n) are the unnonnalized output phase 
samples. 

The recursive normalizer 506 input signals 531 
and 532 and outputs sine and cosine signals 320 and 322. 
The recursive normalizer also outputs recursive feedback 
signals 550 and 552 to complex multiplier 504. The 
recursive normalizer 506 includes two squaring elements 
54 0 and 544, a subtractor 542 , a shifter (divide by 2) 
543, two multipliers 546 and 548, and a constant signal 
"3" 534. The recursive normalizer 506 normalizes 
signals 531 and 532 to very near unity. An exact 
normalization would require dividing the increment phasor 
by its magnitude, which would require a division 
operation, square root, two multiplies, and an add. This 
computation can be replaced by an approximation, which is 
very accurate since the magnitude is already very near 
unity by design. The normalization is defined as the 
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first two terms of a series expansion of the normalization 
factor about one: ^ 



25 



3-J 2 -g 2 
2 



The division by 2 preferably is implemented by shifter 54 3 

Thus, in summary, the output of NCO 314 can be 
expressed as: 

a(n) +j-b(n) - [c (n) +j-d(n) ] -5 

ri^xe cWj-ddrt - Ca( a -i)^*(n-l)l ' [1-0.62^*^1. 
It is important to note that the normalization is done 
within the recursive loop in oscillator 314, and that the 
previous exponential output sample is a factor in the next 
output, which is also normalized. The closed loop in 
oscillator 314 ensures that the magnitude converges to 

unity. . . w ^ 

The following is an analysis of the operation of 

the NCO 314. A bound on the worst case amplitude errors 

can be analyzed by first observing a worst-case increment 

phasor deviation from unity. Setting the derivative of 

the increment phasor magnitude as a function of <P 390 

equal to zero, and then solving for +, will yield the 

maximum increment phasor magnitude error. 



(14) 



d\a(n) +i'b(n)\ s Q . then 4,-0 . 812 

When * is the worst case value, this corresponds to a 
magnitude of 0.951833 (it should be nearly 1). When the 
normalization is applied, the magnitude becomes 0.9965575. 
If this peak error is uniformly distributed about nominal 
unity, then the total mean square error (MSB) is at least 
-60 dB from the signal power. As shown in Fig. 7, the 
actual performance, however, is much better than this 
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bound. Additionally, it is likely that a portion of this 
MSB power is contained at the desired frequency lowering 
the effective MSE. 

As discussed above, setting the increment phasor 
5 to unity is more important than the resulting frequency 

error when setting the real and imaginary elements 518 and 
519. Therefore, a deterministic frequency error exists. 
Although not all NCO applications require absolute 
frequency accuracy, another implementation of the 
10 invention includes an optional portion that corrects the 
frequency. If the desired phase increment is 6 radians 
per sample, then use the following approximation: 

(15) 



/nN 8 ^ u „ a* -1 + v'l +4 -0.6 2 -tan* (0) 

tan (6) = — ; Chen * = g o>t- ar . ta\ \ 

1-0,62 • <& 2 2*(0 . 62xan Id) ) 

In another implementation, the value of <p does not change 
15 nearly as often as the sample rate, so it is practical to 
compensate <p only when needed. An efficient approximation 
for determining a frequency correction (called 
"prewarping" ) is : 

(16) 

<t> n m e a -[l-0.286-eS-(l-0.9-eJ-[l-0.579-eJ]] ] 

20 In an alternate implementation , if needed, arbitrarily 
higher accuracy and/or lower amplitude noise can be 
obtained in NCO 314 by expanding the two series 
approximations to more than the present two terms each. 
One additional imaginary element substantially improves 

25 the accuracy, if the time for additional processing 
exists • 

This technique can directly accommodate a phase 
increment range of at least +-45 degrees (+- 0.125 
normalized frequency). In another embodiment, an entire 
30 frequency range up to the alias can be easily synthesized 
by using phase rotation in quadrants +-(l+-j)» 
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HI. ^.frhud Fil^r Pair Symbol Synchronizer 

Pig. 8 is a schematic block diagram of the 
matched filter pair and symbol synchronizer 240 and symbol 
synchronizer estimator 250 of Fig. 3, which converts a 

5 digitally modulated analog waveform into a digital data 
sequence. The waveform usually consists of filtered 
pulses (e.g., Nyquist filtered for optimal performance in 
noise) whose polarity represents the data stream. The 
elements 240, 250 find an optimal symbol timing for 

L0 sampling an incoming data sequence. The input signals 330 
and 332 are expressed as: 

(17) 

r ( Jc/ f 3 ) «X e i d " %n v ( k/ f 3 -nT) -e* 1 "'^ *n ( kf f a ) 

where d«0, 1, 2, or 3, dependent upon the data pair 
modulating the nth symbol, time t is replaced with sample 

15 times k/fs, p(t) is the matched filter impulse response 
(discussed in more detail below),. *(t) is * random phase 
process (e.g., Rician fading) with an arbitrary mean, fs 
is a sample rate, and n(t) is a complex AWGN process with 
single sided power spectral density No/ 2. 

2 0 The symbol synchronizer 240 includes a pair of 

matched filters 802 and 804. Matched filters 802 and 804 
receive a 3-bit interpolation control signal 872, which 
determines an effective advance or delay of the impulse 
response timing of the matched filters 802 and 804 

25 relative to the input signal 330 and 332. Figs. 17 and 18 
show an input to matched filters 802 and 804. Each 
incoming symbol is represented by, e.g., 8 sequential 
samples. Filters 802 and 804 each contain 8 sets of 
interpolation coefficients for each sample. Thus, the 

30 range of the interpolation control signal is 0 to 7 values 
where this range is modulo one sample period. In the 
described embodiment, a select bits element 870 selects, 
e.g., the 14th through 12th bits of a signal output from 
integrator 860 to be the 3-bit interpolation signal input 
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to the matched filters 802 and 804. Symbol 
synchronization is established when the interpolation 
control sets the matched filter timing coincident with the 
incoming symbol timing. The matched filter timing is 
5 adjusted by using one of e.g., eight sets of filter 
coefficients in each filter , as described below. The 
derivation of the weighting of the filters 802 and 804 
also is discussed below. „. 

Fig. 9 shows more detail of the matched 
10 filters 802 and 804. The matched filter output signals 
803 of matched filters 802 and 804 are: 

(18) 

JC/2 

f**ziy ^ > 3 E < ) / f 3 ) 

Jt— JC/2 
JC/2 

fpunct < n > 3 £ ( ( k+mn) / f s ) 

JC/2 

m *Wi* r < (k+mn-6) / f 3 ) 

JC/2 

where m is the number of complex samples per symbol period 
(e.g., 8), Wkl+i are the filter weights (I is the number 

15 of interpolated filter coefficients between input samples , 
i is the particular interpolated set as indicated by the 
interpolation control signal 872), fs is the sample rate, 
and 5 is the early or late time difference (set to one 
input sample time in the described embodiment) . The 

20 matched filters 802, 804 generally span several or more 

symbol periods to minimize ISI (intersymbol interference) 
while approximating root-Nyquist characteristics. 

The effective sample rate of the matched filter 
relative to the symbol rate determines the resolution of 

25 timing control over the symbol period (i.e., sample points 
per symbol). The more sample points per symbol, the more 
the timing resolution is improved. A key feature of this 
matched filter is its ability to provide very fine 
resolution through interpolation without increasing the 

30 sample rate. This interpolation places virtual sample 

points between the actual sample points. For the filter 
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defined by equation (18), the index n sets the input 
sample index, while the index i sets the interpolation 
index. For example, if there are 8 samples per symbol 
(T-8), a resolution of 64 virtual samples per symbol can 
5 be achieved by using 8 sets of matched filter coefficients 
(i«8), each set skewed by (1/64) of a symbol period. The 
punctual signal is an information bearing signal and is 
passed out of the synchronizer. The early and late 
signals are used to adjust the timing of the filters as 
10 described below. The early and late matched filter 

outputs are rectified through appropriate non-linearity 
elements such as squaring elements 806 in symbol synch 

estimator 250. 

There are three outputs (i-e., early, punctual, 
15 and late) for each I and Q signal path. Hence, as shown 
in Fig. 9, there are actually six matched filters. 
Fig. 17 shows an example of the timing of the matched 
filters leading timing of the incoming symbol, and the 
effect of changing to a different set of the eight sets of 
20 filter coefficients in each of the filters 902 in 

accordance with interpolation control signal 872. Fig. 18 
shows an example of a sample of the timing of the matched 
filters lagging the timing of the incoming symbol, and the 
effect of changing to a different set of the eight sets of 
25 filter coefficients in each filter 902 in accordance with 
interpolation control signal 872. 

As shown in Fig. 8, processor 38 sums the early, 
outputs in adder 810 and sums the late outputs in 
adder 812. Next processor 38 adds the sums in adder 814 
30 and finds the differences between the sums in difference 
element 816. Element 820 performs a normalization by 
dividing the sum by the difference. The output of 
element 820 is an error signal 821 indicating whether the 
matched filter symbol timing is ahead or behind the 
35 incoming symbol timing. 

As shown in Fig. 9, filters 802 provide early, 
punctual, and late outputs for input I. Filters 804 
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provide early, punctual, and late outputs for input Q. 
The three-bit interpolation control signal 872 (also 
called -l" in Fig* 9) selects one of eight sets of filters 
coefficients in each filter 902. The derivation of the 
5 filter coefficients are discussed below. Delay elements 
904 , 906, 904', and 906' each delay the incoming signal 
330, 332 by one sample period. 

Fig, 19 shows the matched filter timing. Each 
of the early, punctual, and late filters 9 02 uses 40 

10 filter coefficients (one of eight possible sets of filter 
coefficients for that filter). If the punctual filter is 
centered on a symbol (solid vertical line), the outputs of 
the early and late filter outputs are equal. If the 
punctual filter is not centered on a symbol (dotted 

15 vertical line), the early and late outputs are unequal. 

Because error signal 821 is a function of the 
square of the input signal level , the loop gain is very 
sensitive to the input level. The first-order gain 
parameter a and the second-order integration parameter fi 

20 set the acquisition time and the tracking bandwidth of the 
closed loop. Additionally, another integrator 860 has 
proved to be beneficial in decreasing acquisition time 
while minimizing tracking (noise) bandwidth. 

The parameters a and R establish an acquisition 

25 time and loop bandwidth which are valid for only one fixed 
input signal level. The effects of these parameters 
change dramatically with varying input levels. The 
dynamic range of the input signal can be extended 
arbitrarily by the normalization described next. 

30 Normalization is achieved by dividing the difference of 
the sum pairs from elements 816 and 818. This 
normalization extends the dynamic range indefinitely while 
maintaining the same acquisition time and tracking 
bandwidth. This assumes the SNR (Signal to Noise Ratio) 

35 is sufficiently high and that finite precision arithmetic 
has no effect* 
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In a preferred embodiment where the signal 
temporarily dips below an unacceptable SNR (e.g., due to 
Hician fading or Rayleigh fading), it is to 
suppress the loop feedback to allow the symbol t^ng to 
"flywheel" through the fades. This fading 
technique is easily implemented by adding a cons tant term 
Eb/N to the normalization denominator. The constant term 
is added at sum element 814. The value of this constant 
term is equal to the sum of the sum pair values when the 
signal level is at the desired break point. The output of 
the normalization in element 820 yields: 

Xin) 3 \f wiy in) \ 2 l^c. in)?+f(E b /N 0 l±mit:) 

where Eb/No is approximately 3 dB over a predetermined 

noise floor. 

The output of integrator 860 is: 



25 



a Jc 

y{n) ia^Vp^ii) 



where gain constant ia a and one of the integration 
constants is B. " i-pl-mentations , a and B 

system para.et.rs that can be set by the user. The second 
integrator constant for integrator 860 is assured to be 
unity. Since th. typical useful dynamic range of the 
incoming signal can be high, and since the **»™**^r 
range is th. number of virtual int.rpolat.on sample, over 
one symbol period, e.g., 64, is typically ™ ch .7"" a 
than the dynamic rang., then it prov.. beneficial to us. 
larger number of bits to handle th. dynamic range of the 
incoming signal, and use some smaller number of bits 
interpolation control. The discrete tim. transfer 
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function of the entire tracking closed loop under high 
signal conditions can be expressed as: 



(21) 



-(z-p) +c*o*z 2 

where c is a constant derived from the timing error signal 
5 after normalization over the piecewise linear region near 
the optimum symbol timing* 

Some plots representative of performance of the 
symbol synchronizer 24 0 of Fig. 3 are shown in Figs 10-13. 
Fig. 10 is a graph of the output of normalizer 820. In 

10* Fig. 10, an initial small bias occurs prior to point 1002 
because the filters are not completely centered on the 
symbol yet. Fig. 11 is a graph of the output of 
integrator 830. In Fig. 11 , symbol synchronization is 
achieved at approximately point 1102. Fig. 12 is a graph 

15 of the output of integrator 840. Fig. 13 is a graph of an 
interpolation control signal 872. In Fig. 13 f 
interpolation control signal 872 dithers over range 1304 
between two output values in an effort to keep the filters 
synchronized with the input signal. Figs. 10-13 relate to 

20 a near worst case example, where the nominal SNR is 6 dB 
(Eb/No), and includes Rician fading (K-10 dB, 200 Hz BW)* 
and other channel impairments. The plots show that even 
under these severe conditions, the symbol tracking 
resulted in a symbol tracking error of less than 1% rms. 

25 Although the previous description focuses on 

symbol tracking, a symbol synchronized in accordance with 
the invention can also be used for symbol timing 
acquisition. It can be shown that the tracking loop will 
acquire the correct timing from any initial symbol timing 

30 offset. The tracking constants a and R can be adjusted to 
yield the best compromise between tracking noise and 
acquisition time. Another implementation of the invention 
dynamically changes these parameters, depending on whether 
the symbol synchronizer is in an acquisition or tracking 

35 mode. In some packet data communications systems, it is 
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customary to provide a preamble clock sequence at the 
beginning of a packet. A typical clock sequence has ones 
and zeros to aid in symbol acquisition. This P* rtiCU ^ 
sequence has a spectrum of mainly line elements on «*tar 
side of the actual carrier frequency spaced at half the 
symbol rate (-1/2T, . Therefore, the signal during this 
clock sequence can be represented by: 

s(n) HA {e :M2 * < '^ /2)flr+ * ) *e Jt " (fe " /l/2,flT ** ) l 
The symbol clock can be produced by forming the 
product of each symbol with its own complex conjugate. A 
phasor containing the phase of this clock can be extracted 
by producing a correlation sample (for a single point m 
time) through the following: ^ 

eJ * = ]£ B~ J2 * tw **\.8in-k) -s' {n+k) ] 

=2A 2 "[i+cos (2nfyir+4>J 
15 The phase of this result can be determined by computing :^ 

.[ «<••* ♦)_" 



$*ARCTAN 



This phase <2» radians = 1 symbol) is added to «-J*-bol 
time modulo one symbol relative to the timing of the input 
sample set. This value sets the initial interpolation 
t Jng control for the matched filter 240- by 
the accumulator of the second integrator with the symbol 
acquisition estimate, effectively receiving instant symbol 
synchronization . 
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IV. Derivation of a Matched Filter 

The following paragraphs discuss derivation of a 
matched filter with a root cosine Nyquist characterization 
with rolloff parameter (b-0% to 100%). 
5 The channel consisting of the cascade of both 

the transmitter and receiver matched filters , and ignoring 
other impairments is to have a raised cosine spectrxjn with 
a roll-off factor of b=60%. The matched filters have a 
spectrum that is the square root of the channel 
10 characteristics. This 60% root-cosine filter is most 
easily represented as follows: 

(25) 



20 



r = 8 b » 0,6 j = y 3 ! 



H(f) = cos 



2-T 2'T 



1-b 1-b 

H(f) =1 ~ ~~2~°T ~2^T 

-os[(i^).(^)] 

The impulse response as a continuous function of t can be 
found by computing the inverse Fourier transform of the 
15 matched filter's transfer function. The impulse response 
is: 

(26) 



I6«jb 2 *t*sin 
hit) »- 



it -(1-b) --|| 4'b'T-coa it-(l+b) --|| sin|*-(l--p: 



tc -(T 2 -16 -b 2 'C 2 ) n-(T 2 -16-b 2 't 2 ) n ' c 

where, e.g., t (in units of sample periods) has the 

following values: 

t - -19 1/2 ... + 19 1/2 for the punctual signal 
t - -18 1/2 ... +20 1/2 for the late signal 
t - -20 1/2 ... + 18 1/2 for the early signal 

Thus, as shown in Fig. 19, the early and late signal lead 

and lag the punctual signal, respectively, by one sample. 
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Other embodiments will be apparent to those 
skilled in the art from consideration of the specification 
and practice of the invention disclosed herein. It is 
intended that the specification and examples be considered 
as exemplary only, with a true scope of the invention 
being indicated by the following claims. 
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What is Claimed i3 : 

1. A demodulator for a mobile 
telecommunications system, including a carrier tracking 
loop receiving a series of symbol samples, comprising: 

a complex multiplier receiving the sequential 

5 received samples; 

a matched filter pair, connected to the complex 
multiplier and outputting a complex filtered signal? 

means for determining a phase difference between 
two sequential received samples; 
10 means for determining a normalized angle in 

accordance with the phase difference; 

means, connected to the normalized angle 
determining means, for determining an increment phasor 
from the normalized angle; and 
15 means, connected to the increment phasor 

determining means and to the complex multiplier, for 
generating a complex exponential value that is input to 
the complex multiplier, 

wherein the complex multiplier, the matched 
20 filter pair, the means for determining a phase difference, 
the means for determining a normalized angle, the means 
for determining an increment phasor, and the generating 
means form a frequency locked loop. 
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2 The demodulator of claim 1, wherein the 
generating means is a numerically controlled oscillator. 

3 A demodulator for a mobile 
telecommunications system, including a ^"^^^^ 

• „ * 9eriea of symbol samples, each having an 
loop receiving a series or sp""* *• 

inphaee signal and a quadrature signal, comprising: 

a comple* multiplier receiving the seguentxal 

received samples; 

a matched filter pair, connected to the complex 

multiolier and outputting a complex filtered signal; 
multiplier^ J & ^ dif ference between 

two sequential received samples; 

means for determining a normalized angle in 
accordance with the phase difference; 

means, connected to the normalized angle 
determining means, for determining an increment phasor 
from the normalized angle; and 

means, connected to the increment phasor 
determining means and to the complex multiplier, for 
generating a complex exponential value that is input 

the complex multiplier, 

wherein the complex multiplier includes: 
means for inputting the complex exponential 
value to the complex multiplier, and 

means for performing a complex multiplication 
operation on the complex exponential value and a current 
pair of the complex inphase /quadrature signals to 
translate the complex inphase /quadrature signals to 
approximately baseband frequency. 

4. The demodulator of claim 3, wherein the 
means for determining an increment phasor includes 
normalization means for allowing the carrier tracing loop 
of the demodulator to flywheel through a signal fade. 

5. The demodulator of claim 4, wherein the 
normalizer includes a NORM constant. 
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